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I , Abstract. We review statistical theories and numerical methods employed to 

consider the sample size dependence of the failure strength distribution of disordered 
materials. We first overview the analytical predictions of extreme value statistics and 
\ fiber bundle models and discuss their limitations. Next, we review energetic and 

geometric approaches to fracture size effects for specimens with a flaw. Finally, we 
overview the numerical simulations of lattice models and compare with theoretical 
J> ' models. 
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1. Introduction 

Understanding how materials break is a fundamental problem of science and engineering. 
The difficulties stem from the non-trivial dependence of the fracture strength on the 
characteristic lengthscales of the samples. This was already noted by Leonardo da 
Vinci, who measured the carrying-capacity of metal wires of varying length pQ. He 
observed that the longer the wire, the less weight it could sustain. The reason for 
this behavior is rooted in the disorder present in the material. If a sample can be 
divided into non-interacting subvolumes, then the strength is dominated by the weakest 
one and its distribution could in principle be computed using extreme value statistics 
[2J, [3] . Longer wires are more likely to possess weak parts and are thus bound to fail at 
smaller loads on average. The quantitative understanding of this statistical size effect is 
difficult and the low tails of the strength probability distribution are not easy to sample 
experimentally. In addition, microcracks interact by long-range stress fields. Hence, 
the assumption made to use extreme value statistics, that of independent subvolumes, 
is difficult to justify in practice. These problems are particularly compelling in quasi- 
brittle materials, such as concrete and many other composites, where sample failure is 
preceded by significant damage accumulation pE]. 

The typical materials science question is to study the size effect with a pre-existing 
flaw, a notch. Failure in quasi-brittle materials is then determined by the competition 
between deterministic effects, due to the stress enhancement created by notch, and the 
damage accumulation around the defect due to the stress concentration [5] [6] . The effect 
of disorder is then treated in an effective medium sense by defining a Fracture Process 
Zone (FPZ) around the crack tip. Starting from these observations, several theoretical 
formulations based on linear elastic fracture mechanics (LEFM) have been proposed in 
the literature and compared with experiments [5j [6] [7J [8] [9j [10], [11] . The problem is how 
to extend LEFM when disorder is strong and can not be treated as a small perturbation 
on the homogeneous solution. 

The statistical physics approach to fracture is based on simple lattice models, which 
allow for a relatively simple description of disorder and elasticity. The models are 
sometimes amenable to analytical solutions, and usually are simulated numerically. In 
the simplest approximation, elastic interactions are replaced by a load transfer rule which 
is applied when elements fail or get damage. These "fiber bundle models" can be solved 
exactly in some cases and can thus provide a useful guidance for the simulations of more 
realistic models in which the elastic medium is represented by a network of springs or 
beams. In this case, the local displacements can then be found by standard method for 
solving coupled linear equations. Disorder is modeled for instance by imposing random 
failure thresholds on the springs or by removing a fraction of the links. The lattice 
is loaded imposing appropriate boundary conditions and the fracture process can be 
followed step by step, in a series of quasi-equilibria. The cornerstone in this respect has 
been for the last twenty years the Random Fuse Model (RFM) |1 2J , a lattice model of 
the fracture of solid materials in which as a further key simplification vectorial elasticity 
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has been substituted with a scalar field. 

In this paper, we review the statistical physics approaches used to understand size 
effects in fracture. We first discuss the traditional weakest-link approaches and the 
solution of fiber bundle models. Next, we consider size effects arising from energetic 
considerations when a single dominating crack is present, and when its growth may be 
of importance. Finally, we discuss the results of lattice models for fracture and compare 
the results obtained with the theoretical arguments discussed previously. 

2. Size effect from extreme value statistics 

2.1. The weakest link model 

The key theoretical concepts needed to understand the distribution of fracture strength 
and the associated size effects in heterogeneous media date back to the pioneering work 
of Gumbel [2], Pierce [13] and Weibull [3] on the statistics of extremes. The general 
idea stems from a weakest link argument: the failure strength of an extended object is 
ruled by its weakest local subvolume. For a disordered system, the larger the sample 
the easiest it is to find a weak region. 

The easiest way to illustrate this concept is based on a one-dimensional model. 
Consider a chain composed by N elastic links that can sustain at most a stress Sj without 
breaking. We assume that the threshold stresses s« are independent random variables 
distributed according to a probability density function p(si) and a corresponding 
cumulative distribution V{s,i < a) = P{o) = Jq p(x)dx. The chain will fail as soon 
as one of the links, the weakest one, fails. Hence, the probability that a chain fails at 
a stress larger than a is equal to the probability that the minimum of Sj is larger than 
a, which is just the probability that all the links have thresholds larger than a. This 
condition can be expressed in mathematical terms for a chain of N links as 

i - Psifr) = n ><r) = (i- p(°)) N , (i) 
i=i 

where Pat (a) is the probability that the failure stress is smaller than a. In the large N 
limit, Eq. ([I]) can be approximated as 

P N (a) ~ 1 - exp[-NP(a)]. (2) 

Under some broad assumptions on the tail of P(cr), it can be shown that Pn{ct) converges 
to a stable asymptotic form (i.e. a distribution with a shape that does not depend on 
N). In more technical terms a stable distribution is invariant under Eq. (OQ) provided 
the variable is linearly transformed according to cr N = a N a + b N . For large N, the 
statistics is dominated by the low value tail of the distribution P(a). If we assume it to 
scale close to o = as P(<r) — (cr/cr ) M , where a is a characteristic stress. Under this 
assumption, we obtain 

P N (a)~l-e^[-N(a/a )% (3) 
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which is invariant under Eq. §I§ if we transform the variable a with = N 1 ^ and 
b N = 0. 

Eq. ([3]) is the celebrated Weibull distribution [3], predicting that the average 
strength of a chain composed by N links decreases as 

(a N ) oc AT 1 /". (4) 

The Weibull distribution represents still today the main tool used to analyze failure 
statistics in various materials, although the validity of its underlying assumptions is 
in general difficult to demonstrate. Real samples can not generally be schematized 
as a chain of independent elements with random failure thresholds. In many cases, 
such as in quasibrittle materials, the sample does not even fracture at once but 
sustains a considerable amount of damage before failure. Furthermore, long-range elastic 
interactions could correlate different regions of the sample invalidating the assumptions 
used to derive the Weibull statistics. That the interactions and microcrack grow would 
be irrelevant at sufficiently large scale, leading to the recovery of the Weibull distribution, 
has never been proven rigorously. 



2.2. The largest crack model 

It would be desirable to relate the failure statistics to some geometrical characteristics 
of the microstructure of a material, going beyond a simple description of a sample as a 
collection of regions with different random strengths. An attempt in this direction was 
proposed by Freudenthal [H] and then rederived later by various authors for different 
models jT5l [T6l [TTj . The idea is to schematize the microstructure as a set of independent 
linear cracks of length cij, distributed according to given probability density function 
p(aj). If the cracks are sufficiently diluted, we can treat them as completely isolated 
and assess their stability by the simple energetic argument due to Griffith. The idea 
is establish the conditions under which crack growth becomes energetically favorable. 
Considering for simplicity a two dimensional geometry, the elastic energy released by a 
crack of length a is given by 

E* = (5) 
2E 

where a is the applied stress and E is the Young modulus. Forming a crack involves the 
creation of crack surface with an energy cost which in the ideal case would be related 
the rupture of atomic bonds, 

£ SUTf = 2aG, (6) 

where G is the fracture toughness. The crack becomes unstable when the total energy 
decreases 

*L = -7ra 2 a IE + 2G < 0, (7) 
da 

which implies that a crack of length a becomes unstable when the stress is equal to 



2E G 
V vra 
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If we have a collection of N independent random cracks in a volume V, we can 
assume that the sample will fail when least stable one will become unstable. The 
problem thus reduces to find the distribution of the largest crack. As discussed in 
Sec. I2.1[ the cumulative distribution is given by 

P N (a)~exp-\pVP(a)], (9) 

where p = N/V is the crack density and P{a) is the cumulative distribution associated 
with the crack density distribution p(a). The asymptotic behavior of Eq. ([9]) is ruled 
by the large- value tail of P(a). It can be shown that if the tail decays faster than an 
exponential, Eq. (Q converges to the Gumbel distribution [2J for large N 

Pn{u) — exp[— pV exp(— a/a c )], (10) 

where a c is the characteristic scale of the crack length distribution. Combining Eqs. ffTUl) 
and (jSJ), one can derive the strength distribution as 

-pV exp 
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a 

where <7o = 2EG/na c . Eq. (TTTT) predicts a size effect for the average strength of a type 
logarithmic in size: 

(ov)~A/(B + Clog(V)). (12) 
2.3. Interacting cracks: fiber bundle models 

The main shortcoming of extreme value statistics is the fact that local failure stresses 
are considered independent. This is difficult to justify in practice since cracks induce 
long-range interactions that may correlate stresses in different regions of the sample. 
A simple way to analyze the role of interaction in the fracture of disordered media is 
represented by the study of fiber bundle models [131 EH]- These models consider a set 
of brittle fibers, with random failure thresholds, loaded in parallel. When a fiber breaks 
its load is redistributed to the other fibers according to a prescribed rule. The simplest 
possibility is the case of an equal load sharing (ELS), in which each intact fiber carries 
the same fraction of the load. This case represents a sort of mean-field approximation 
and allows for a complete analytic treatment [T8l [19j [20j [21], [22] . At the other extreme 
lies the local load sharing model (LLS) where the load of a failed fiber is redistributed 
to the intact neighboring fibers [231[2l[25l[2a[2?l[^ 

We consider first the case of ELS fiber bundles, in which N fibers of unitary Young 
modulus E = 1 are subject to an uniaxial load F. Each fiber % obeys linear elasticity 
up to a critical load X{, which is randomly distributed according to a distribution p(x). 
When the load on a fiber exceeds Xi, the fiber is removed. Due to the ELS rule, when 
n fibers are present each of them carries a load = Fjn and consequently a strain 
e = Fjn. The constitutive law for ELS fiber bundles can be easily be obtained from a 
self-consistent argument. At a given load F, the number of intact fibers is given by 

/ r F/n \ 

n = N 1 - / p{x)dx . (13) 
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Rewriting Eq. (1131) as a function of the strain, we obtain the constitutive law 

l = e(l-P(e)), (14) 

where P(x) is the cumulative distribution obtained from p(x). For simplicity we rewrite 
Eq. as / = e(l — P(e)) where // = F/N is the load per fiber. Failure corresponds 
to the maximum e c of the right-hand side, after that there is no solution for e(f). 
Expanding close to the maximum we obtain / ~ f c + B(e — e c ) 2 , which implies that the 
average rate of bond failures increases very rapidly before fracture: 

%~U>-!T». ("0 

Contrary to the case of extreme type models, the strength of a fiber bundle with ELS 
does not vanish in the large N limit. In particular, for any threshold distribution such 
that 1 — p{x) goes to zero faster than 1/x for x — > oo, the strength distribution is 
Gaussian with average f c = e c (l — p(e c )) and standard deviation <jf = e c p(e c )(l — 

p(e c ))/VN m- 

A typical LLS model considers a one dimensional series of fibers loaded in parallel 
with random breaking thresholds from a distribution p(x). When the load on a fiber 
exceeds the threshold its load is redistributed to the neighboring intact fibers. Thus 
the load on a fiber is given by /j = /(l + k/2), where k is the number of failed 
fibers that are nearest neighbors of the fiber i and / = F/N is the external load [23J. 
Even for this apparently simple one dimensional model a closed form solution is not 
available, but several results are known from numerical simulations, exact enumeration 
and approximate analytical methods [23l [2H ESI [26l [271 12HI [291 EOl EH [32l [33] Contrary 
to the ELS model, LLS fiber bundles normally exhibit non-trivial size effects as could 
be anticipated from general consideration of extreme value statistics. In particular, in 
the limit of large N, it has been shown that the strength distribution follows the form 

W(f) = 1 - [1 - C(f)] N ~ 1 - exp(-NC(f)), (16) 

where C(f) is a characteristic function, close to the Weibull form, but difficult to 
determine exactly [23J. The existence of a limit distribution has been recently proved 
under very generic conditions for the disorder distribution in Ref . [31] , but a numerical 
estimate of its form typically requires extremely large samples. From Eq. [16] follows 
that the average bundle strength decreases with N as 

/ c ~l/log(JV), (17) 
so that an infinitely large bundle would have zero strength. 



3. Size effects from energy and geometry 

In the previous section, we discussed the statistical size effect, which is caused by the 
randomness in the material strength. In this section, we discuss two alternate ideas 
relevant, namely, the energetic size effect and the geometric size effect on the failure 
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strength of quasi-brittle materials. The energetic size effect is based on the Griffith's 
criterion of energy balance for stable crack propagation. Intuitively, the energetic size 
effect arises due to stress redistribution and the associated stored energy release as a 
large fracture process zone (FPZ) develops ahead of the crack tips. The size of the FPZ 
introduces an additional length scale into the problem and its influence on the fracture 
strength is the energetic size effect in quasi-brittle materials. This idea does not derive 
from statistical fracture, but turns out to be quite useful in understanding the results 
of the next section. 

On the other hand, the geometric size effect is a consequence of fractal geometry 
of crack surfaces. The morphology of crack surfaces in disordered media can then be 
described by a self-affine transformation, which implies a length-scale dependence of the 
surface energy of a crack. This effect is supposed to introduce a geometric size effect on 
the fracture strength. 

3.1. Energetic size effect 

Classical continuum theories of elasticity, plasticity and damage mechanics do not 
possess a characteristic material length scale I. Consequently, the nominal strength 
of geometrically similar structures is independent structure size, i.e., there is no size 
effect. Alternatively, the presence of a characteristic length scale introduces a size effect 
on the nominal strength of the material. 

Recall the basic idea of Griffith's criterion, that during quasi-static crack 
propagation the available energy must be equal to the energy required to create new 
crack surfaces. This is precisely the expression derived in Eq. (jSJ) for brittle materials 
under plane stress criterion. For quasi-brittle materials, i.e., when the fracture process 
zone is comparable to or larger than the structural size, Griffith's criterion can be 
generalized using an effective crack size concept that accounts for the influence of FPZ 
size on nominal strength. 

Following the classical continuum fracture mechanics treatment, a generic 
expression for the stress intensity factor around a crack of size a may be expressed 
as 



where k{q) denotes a dimensionless function, a = a/L denotes the relative crack size 
with a = ao + Aa such that ao is the initial crack size and Aa is the incremental crack 
size, and L is a characteristic system size. Equation (fT8j) represents an extension of 
the classical stress intensity factor around a central notch of size 2a in an infinite panel 
(Ki = a^fiui with K,(a) = I) to finite system sizes and various boundary conditions. 
The influence of characteristic system size L on stress intensity factor can be seen clearly 
by rewriting Eq. ( fT8i) as 






(19) 



where 0(a) = y/naK,(a). 
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Using Irwin's relation, the generic expression for elastic energy release rate Q can 
be written as 

Q = § = ^9(a) (20) 
where g(a) = (p 2 (a). For quasi-static crack growth, Griffith's criterion states that 

g = n (21) 

where 1Z represents the crack resistance curve (R-curve: TZ(Aa)), which is expressed 
more generically as 

K = R ip(a). (22) 

As mentioned earlier, the fracture of quasi-brittle materials is preceded by the 

development of fracture process zone, wherein the material undergoes progressive 

damage due to microcracking and interface failure. The presence of the FPZ introduces 

a characteristic length scale I into the problem in addition to the already existing length 

scales: ao, Aa and L. A simple dimensional analysis dictates that the elastic energy 

release rate and crack resistance curves be expressed in terms of dimensionless ratios 

Aa , „ i 

T) = — and i? = - (23) 

in addition to the initial relative crack size a = 9 f;. Since a = a + 77$, the energy 

release rate and the crack resistance can be expressed as 
2 7- 

g =^g(a ,r),#) (24) 
K = RTP{a ,r],#). (25) 
For fracture to occur at the peak load, we require that 

OAa dAa 1 ' 

Using the relations G = R and r\ — ^ , Eq. (126]) may alternatively be expressed as 

d(fi(a ,ri,-d) _ dip{a ,r),ti) 
dr) dr) 

From Eq. (j27p . one can derive an implicit expression for the critical crack as 

77 = (28) 

Substituting Eq. ( |28l) into Griffith's fracture criterion g = 1Z and simplifying the result, 
we get 

(29) 

where 

h{a ,$) = — =—r. 30 

^(a ,?7 c (ao,^),^) 



(27) 
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For very large system sizes L, i? approaches zero asymptotically, i.e., $ — > 0. Hence, 
/i(a ,tf) can be approximated as 

fc(a o> 0)«/io(ao)+M (31) 

where 

&o = /*(ao,0) (32) 
_ dfr(a ,?9) 

d$ '* =0 " ^ ' 

Substituting Eq. <^B} into Eq. ([23]), we have 



^ (34) 



VLho + Ihi 

For quasi-brittle materials with an initial notch and constant R-curve, we have 
ho = g(a ) ~ a = a /L. Denoting £ = £h\, Eq. [M]can be specialized to the initially 
notched samples as 

°° = V^T? (35) 

Note that this argument states that £ depends on the partial derivative of the ratio of 
Q and 1Z on so it is related on the relative scales of the two rates. 



3.2. Geometric size effect 

In this section, we discuss the implication of roughness of crack surfaces on the linear 
elastic fracture mechanics and therewith on the scaling of strength. For many years 
now, it is hoped that there exists a simple relation between material toughness and the 
self-affine exponent of rough cracks, although such a hope appears too optimistic now 
(see Ref. [35] for controversies related to this topic). An exception is the indirect fact 
that explaining crack roughness via the theory of depinning of elastic manifolds implies 
that the critical stress intensity factor can be expressed via a few relevant parameters 
such as crack front elasticity and the strength of the disorder, while the precise value 
depends on the geometry (planar cracks, notches or cracks in two or three dimensions...) 
|36[ [3TI [38] 39J. There is, however, no size effect except for the one that arises from 
the finite crack length and the resulting correction to the critical stress intensity factor, 
in analogy with finite size corrections observed in critical phenomena [371 HQ]. It is 
important to note that the roughness and strength can be directly coupled, since e.g. 
the crack growth resistance and the roughness exponent (fractal dimension) follow from 
the same process. In what follows, the empirical observations of strength are explained 
a posteriori using as input a given roughness exponent for the crack. 

In particular, roughness of crack surfaces influences both (i) the stress concentration 
around the crack [H] [42], |43] and (ii) the energy required to create crack opening through 
its dependence on the actual crack area [44]. Experiments on several materials under 
different loading conditions suggest that the crack surfaces exhibit self-affine scaling, 
which implies that if the in plane length scales of a fracture surface are scaled by a 
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factor A then the out of plane length scales (height) of the fracture surface scales by A*% 
where ( is the roughness exponent. This also implies that the actual (curvilinear) length 
of the crack is a non-trivial function of length-scale. Consequently, an infinitesimal 
crack opening will now cost more surface energy than in the usual case, and thus the 
Griffith's criterion is modified in a way to account for the increased energy requirement 
to propagate cracks, and possibly to account for the modified stress fields and hence the 
elastic energy due to a rough crack. 

Following Refs. [HI H2J H3], the asymptotic near-tip stress fields around a rough 
crack may be expressed as 

(Tij = Krr-ffifi), (36) 

where r is the radial distance from the crack tip, f(8) denotes the functional dependence 
of stress concentration around crack tip, and (3 = 1 — l/2( for 0.5 < ( < 1 and (3 = 
otherwise. Equation ( 1361) implies that the stress concentration around a rough crack 
can be expressed as [43] 

Kj ~ aa p X {P) (37) 

where xW) is a dimensionless function. The influence of characteristic system size L on 
stress intensity factor is incorporated into Eq. fl37|) as 

Kj = <j)(a) (38) 



a 



L 2f3 



G = ^-g{a ^$). (39) 

From the crack resistance point of view, it is expected that the roughness of a 
self-affine crack leads to R-curve behavior. Specifically, for a self-affme crack, the crack 
surface energy scales as 

£ surf = 2 7 a 1 ^ (40) 

where 7 denotes the specific energy per unit of fractal measure and 1/C denotes the 
fractal dimension of the crack. Hence, the crack resistance is given by 

R = 2 7 a 1/c ~ 1 (41) 

We also note that a size effect law based on anomalous scaling of crack surfaces has 
been proposed in Refs. [9j [10] . In the case of anomalous scaling of crack surfaces (see 
Eq. ( l4"2l) below), the extra scale also enters the scaling argument, further changing the 
amount of surface energy needed for infinitesimal crack advancement. This implies a 
non-trivial "R-curve" , or crack resistance, since the energy consumed via forming new 
crack surface depends not only on the full set of exponents but again on the distance 
propagated. An expression can be written down using solely geometric arguments and 
the exponents, such as 



Ah(£,x) ~ A 



floc^i-Cloc jf £<£( x ) 

if 
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where Ah(£, x) denotes the height fluctuations of fracture surfaces estimated over a 
window of size £ along the x-axis and at a distance x from the initial notch. Also, 

= Bx x l z is a crossover length along the x-axis below which the fracture surface 
is self-affine with local roughness exponent (i oc . This crossover length depends on the 
distance x to the initial notch. For empirical reasons, one needs to note that there 
should be also a lower cut-off above which the self-affinity can be observed |45j . and 
that the argumentation should change if the crack geometry is very branched [4"S] . 

There are two different regimes: for £ 3> £(x), Ah(l, x) ~ x^, and for i <C it 
is characterized by the exponent (( — (i oc )/z, where z is the dynamic exponent. Since 
the effective area is (from which surface energy follows) 



where 8A r is the real area increment and SA p its projection on the fracture mean plane. 
Morel et al. make now the simplest first order approximation, that Eq. ( 1421) can be 
approximated as G = 2 , ys(x)/L, where s(x) is the length of the crack profile that can 
be estimated by covering the profile path with segments of length 5 whose horizontal 
projection on x-axis is £q. The elastic energy released is not taken to be affected by 
the self-affinity of the crack. Consequently, we have s(x) = (L/£ )5, which can be 
approximated as 



Substituting Eq. H3]into G = 2js(x)/L, we get an expression for crack resistance 

as 



within the region of crack growth (Aa <C Aa sat = x sat ), and when Aa ^> Aa sat , the 
crack growth resistance becomes 



4. Size effects from disordered fracture models 

Simulations of statistical fracture provide a tool to analyze the validity of various 
scenarios for the size-dependence of strength. The major tool for this has been the 
Random Fuse Model [121110], where one approximates continuum Linear Elastic Fracture 
Mechanics (LEFM) with a spatial discretization (lattice) and a scalar analogy (voltages) 
for the displacements. The material response is formulated in terms of "fuse elements" 
which have a conductivity Go and a failure threshold i c . Usually the "elasticity" or 
the fuse conductances is kept constant, while various threshold distributions p(i c ) are 



G 5A P = 2 7 8A r 



(42) 




x y§> X sa t 
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employed. Until now almost all the effort has been spent on two-dimensional systems 
for the sake of numerical convenience. The discussion of established results can be 
divided into two topics, studies of the weakest-link effect and analyses of the size effect 
for notched specimens. 

4-1. Statistical size effect 

The major results on the strength of the RFM were obtained in the late 80 's by Duxbury 
and co-workers [I7J, H7J . The starting point is to consider weak disorder employing a 
failure threshold distribution composed by two delta-functions at i c = and i c = 1 
which naturally translates into " dilution disorder" or the removal of a fraction 1 — p of 
the fuses at the beginning. In the diluted limit (p ~ 1), it is possible to formulate a non- 
rigorous theory of the size-effect. First, the flaw size statistics indicates an exponentially 
decaying length distribution for the microcracks that are induced into the samples by 
the dilution disorder. Second, in practice the critical current J c essentially corresponds 
to the current I\ st which the weakest fuse fails. This means that the case under study 
is very brittle, and the crack growth resistance is small. Given these two assumptions, 
Duxbury et al. write down a modified extremal statistics theory of the type of Eq. 
( fTOl) . and compare it to RFM simulations in two dimensions with small dilution. There 
are straightforward generalizations of the weak-disorder dilution theory towards p ~ p c , 
the bond percolation of the RFM geometry (for square lattices, p c = 1/2), since at the 
percolation point the outcome can be simply understood via the "red bond" or singly- 
connected bond arguments of percolation theory [37]. This limit is not relevant for most 
materials, but in general the defect size statistics can not be assumed exponential, and 
often a noticeable damage accumulation takes place before failure. 

To this end, we consider a variant of the RFM where p(i c ) has a finite support which 
extends down to zero. In particular, we consider disorder with a cumulative distribution 
P(ic) — V) z c G [0,1]. The case A = 1 corresponds to the uniform distribution, 
while different values of A allow to tune the strength of disorder. Since a fuse can fail 
even at low currents, with these kind of distributions there is a considerable damage 
accumulation before failure. We illustrate this in Fig. [T] which shows the relation 
between the accumulated damage d and the failure stress a c = I c /L The figure shows 
the results of two dimensional simulations for three different disorder strengths A and for 
several system sizes from L = 64 to L = 320 or even 512 (for the smallest A or weakest 
disorder). It is interesting to note that the d(a c ) relation follows roughly a power-law. 
While this general behavior could be anticipated on the basis of a simple mean-field 
relation d = J CTc di c p{i c ) = cr^, the measured curves do not follow this prediction. As 
we discuss below in more detail, the damage is a sum of a statistically homogeneous 
background plus the FPZ contribution so the failure of the mean-field is not a surprise. 

Damage accumulation is responsible for a fundamental complication in the 
evaluation of size effects which should be relevant for real materials. To relate this to the 
discussion in Section 2, the cumulative P(a) in Eq. (Q becomes dependent on a. Fig. [2] 




shows numerical results for the cumulative strength distribution P(<J C , L) for the RFM. 
As demonstrated in the Figure, the data are not described by a Weibull distribution, 
but not even the modified Gumbel distribution mentioned above captures the observed 
behavior. It is interesting to note nevertheless, that for larger L the double exponential 
plot seems to reveal the presence of an asymptotic scaling form. More importantly, 
the behavior here shows that the results are not easily interpreted in terms of extreme 
statistics in its simples forms. It would be desirable to formulate a theory capturing the 
evolution of P(i c , L) as L increases. Studies of idealized chains-of-fiber-bundle models 
have indicated the presence of cross-over effects depending on system size and disorder 
and the failure of standard extremal statistics scalings [HJ H9[ [50]. Such chains consist 
of simple one- or two-dimensional models in series. Thus their behavior is dictated 
when V changes, in analogy to Sec. 12.11 by the chain length N and the distribution 
-P(cr). Important physics originates at the low-strength tails of the element (fiber, fuse, 
bond, element) strength distribution and from the mesoscopic damage dynamics. Other 
important effects are crack arrest by strong elements, and for weaker disorder the concept 
of a critical defect size, reaching which anywhere in the sample leads to catastrophic 
failure. 

4-2. Strength of materials with flaws 

Another fundamental question is whether Eq. (jSJ) from the Griffiths' argument remains 
valid for realistic materials. A suitable generalization as discussed in the previous section 
is a c = K c j\J^ + do, where we have added a scale £ and written K c ~ y/EG c for the 
fracture toughness (numerical factors relating eg. plane stress or plane strain loading 
scenarios have been omitted). £ plays here the role of incorporating a disorder effect: 
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Figure 2. Cumulative strength distributions for the RFM (with A = 0.6) plotted for 
several L. Upper panel: the data shown using a "Weibull paper". If the data would 
obey the Weibull distribution, they should all collapse into a single a straight line. 
Lower panel: a generalized "Gumbel" plot. Again a collapse to a straight line would 
be required if the data were described by a generalized Gumbel distribution. 



that cracks are masked when they become very small by fluctuations, and thus the 
strength saturates to K c j\f^. As noted by Bazant (see e.g. [8]) a natural reason for the 
scale-length is the Fracture Process Zone which is non-negligible in quasi-brittle media. 

Extensive numerical simulations of the RFM and also other models - Random 
Spring Model, Random Beam Model - have allowed to understand the physics of size 
effect here and to elaborate on a scaling theory [TTI [51] . Figure [3] shows an example 
of 2d RFM data plotted using the inverse square of Eq. f)35p . For small flaws one 
sees a cross-over away from the behavior predicted by the above scaling, and the inset 
illustrates how this depends on the specimen size, which is otherwise absent in principle 
from Eq. ( 1351) . 
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Figure 3. RFM strength data for intermediate disorder and varying flaw sizes ao and 
sample linear sizes L. The inset compares fracture strength in unnotched samples with 
that predicted by Eq. ([55]) for ao = 0. 



The dependence of the size-effect on L and the disorder strength A (in full 
generality, we take A = to correspond to pure systems with no disorder) can be 
summarized with the scaling 

= £ + a o/(a c /a ) (43) 

® c 

where the scaling function f(y) follows 

f(y)-\ . f 44 
1 y if y » i 

The important point here is that Eq. (I44p takes into account of the weak-link effects 
discussed above. The length-scale a c ~ (K c (A)/a(L, A)) 2 — £(A) indicates the cross-over 
between that regime and the LEFM one where Eq. ( |35l) is valid. 

The size-effect can thus be summarized with the aid of two parameters, £ and 
K c = \J GE where E is the elastic modulus, and G the fracture energy. It has to be 
emphasized that all the three quantities of £, E, and G depend on the presence of 
disorder. Figure H] illustrates the FPZ damage profile, which underlies £. It measures 
the exponential decay of the damage from the crack tip, in a situation where the LEFM 
stress profile (and its correspondence in the RFM) is not seen in the damage due to 
screening. The stronger is the disorder, the larger is £ as one would expect, indicating 
also a strength reduction. It would be very interesting to develop an analytical model of 
this screening effect. The figure H] demonstrates however that there are strong sample-to- 
sample fluctuations: the exponential damage profile and its role in the average strength 
behavior or size-effect is to be understood as a statistical average. 

It is interesting to note that these statistical models exhibit a non-trivial R-curve 
if one computes the crack extension using the changing scale of the damage cloud, £ (a) . 
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Figure 4. The damage at maximum stress a c for a RFM with linear size L = 128 
and intermediate disorder A = 0.6. The black markers illustrate the broken fuses in a 
single sample, and the color code depicts the average damage over N — 2000 samples. 
One can see screening effects next to the notch boundaries (weak damage), and then 
two clouds of damage next to the notch tips. 



The strength prediction of Eq. fl35l) is then "self-consistent" . For those ao for which it 
is valid, £ is a constant to which value £(cr) grows as damage accumulates from zero. 
Likewise, K c is to be measured at a c as well. There are no clear signs of any kind of 
geometrical size effects, since though the cracks formed are self-affine, the growth takes 
place after the stress maximum [51] . However, the damage cloud is truly two-dimensional 
and we thus note that the results differ from those obtained for uniaxial fiber composites 
[521 [53] which exhibit a scaling similar to Eq. ( 1351) but with a logarithmic correction and 
which reproduce a R-curve in a similar fashion as the statistical models. 

This brings us to the important issue of strength fluctuations in the presence of 
a dominating notch. In the previous sections we have discussed the important role 
of extremal statistics in this. Now, the size-effect, of the whole distribution, is not 
expressed as a function of L but of a . Figure [5] illustrates the behavior of P(a c ,a ). 
As ao grows the average strength goes down, and the width of the distribution reaches 
a maximum at an intermediate value of ao- For small notches, ao < a c the distribution 
naturally is close to the finite-size form for p(a c , ao = 0, L). For large values of ao we 
see a narrowing of the distribution (in fact, Aa c /a c goes to zero with a increasing) and 
the distribution is well approximated in the central part by a Gaussian shape. This 
is interesting, since the Gaussian distribution is also obtained in global load sharing 
fiber bundles [18] and in chains-of-fiber-bundle models with local load sharing and 
heterogeneous fiber strengths [50J. However, even for the largest ao = 48 the tails 
are broader than in the Gaussian case. It would perhaps be expected that for the weak 
strength limit the fuse threshold distribution would be relevant [H], and for this reason 
it would be interesting to study these models also with e.g. a Weibull P(i c )- 
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Figure 5. The distribution p(a c ,ao) for a RFM with linear size L — 192 and 
intermediate disorder A = 0.6. Three different a are used: L = 4,16,48. The 
last dataset is compared to a Gaussian fit in the central part around a c — (cr c ). The 
numerical distribution are obtained from 1000 to 4000 samples depending on A and 



5. Conclusions 

In this paper we have reviewed various statistical approaches to fracture size effect. We 
first discussed the Weibull theory and similar approaches based on the weakest link 
concept. According to these extreme value statistics arguments, larger samples should 
be more likely to possess weak regions and are thus bound to fail at a smaller load. While 
this general statistical argument is generally believed to be correct, the presence of the 
stress enhancement around a crack and damage accumulation prior to failure make a 
quantitative theory for size effects in materials still elusive. In cases when a dominating 
crack, or a notch, is present one can rely on energetic arguments that link the crack 
geometry to the failure stress. Further complications arise when the dominating crack 
is not straight but has a self-affine geometry, as observed in experiments. 

These general statistical approaches have been tested using statistical lattice models 
for fracture, where a set of discrete elements with random failure thresholds are subjected 
to an increasing load. These models allow to study the crossover between notch 
dominated size-effects and statistically induced ones. Despite this, the role of damage 
accumulation for fracture size effects in unnotched samples still remains unclear. On the 
experimental side, large statistical sampling would be needed to reach firm conclusions 
about the asymptotic behavior of materials strength. The above observations apply for 
quasi-brittle fracture first and foremost, and if viscoelastic or viscoplastic deformation is 
important our understanding is much less developed. These important questions remain 
open for future investigation, so that statistical fracture promise to pose interesting 
challenges for the years to come. 
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